The Paradox of the expected Time until the Next 

Earthquake 



D. Sornette and L. Knopoff 

Institute of Geophysics and Planetary Physics 
University of California, Los Angeles 



1 



Abstract 

We show analytically that the answer to the question, "The longer it has been since 
the last earthquake, the longer the expected time till the next?" depends crucially 
on the statistics of the fluctuations in the interval times between earthquakes. The 
periodic, uniform, semigaussian, Rayleigh and truncated statistical distributions of 
interval times, as well as the Weibull distributions with exponent greater than 1, 
all have decreasing expected time to the next earthquake with increasing time since 
the last one, for long times since the last earthquake; the log-normal and power- 
law distributions and the Weibull distributions with exponents smaller than 1, have 
increasing times to the next earthquake as the elapsed time since the last increases, 
for long elapsed times. There is an identifiable crossover between these models, which 
is gauged by the rate of fall-off of the long-term tail of the distribution in comparison 
with an exponential fall-off. The response to the question for short elapsed times 
is also evaluated. The lognormal and power-law distributions give one response for 
short elapsed times and the opposite for long elapsed times. Even the sampling of 
a finite number of intervals from the Poisson distribution will lead to an increasing 
estimate of time to the next earthquake for increasing elapsed time since the last one. 
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Introduction 



While small earthquakes after removal of aftershocks have a Poissonian distri- 
bution (Gardner and Knopoff, 1974), intermediate and large earthquakes in a given 
region are clustered in time (Kagan and Knopoff, 1976; Lee and Brillinger, 1979; 
Vere- Jones and Ozaki, 1982; Grant and Sieh, 1994; Kagan and Jackson, 1991, 1994; 
Kagan, 1983; Knopoff, et al., 1996); 'clustering' is taken to mean that the earthquakes 
do not have a purely Poissonian, memoryless distribution of time intervals. According 
to this interpretation, periodic earthquakes are the extreme limit of clustering and 
can be predicted exactly. More generally, if there is temporal clustering, the estimate 
of the probability of occurrence of a future earthquake in a given time interval is im- 
proved if there is a knowledge of the times of previous events, since clustering implies 
a memory. To express this property quantitatively, we relate the elapsed time since 
the last earthquake in a region to the conditional probability of occurrence of the 
next earthquake within a given time interval from the present. Davis, et al. (1989) 
have posed a version of this problem in the form of the question, 

(Q.) (can it be that) "The longer it has been since the last earthquake, the longer 
the expected time till the next?" 

The observation of Davis, et al. for the log-normal distribution was that the answer to 
Q. is positive. Ward and Goes (1993) and Goes and Ward (1994) showed numerically 
that, in the case of the Weibull distribution, the response to Q. can be either yes 
or no, depending on the exponent in the distribution. The positive responses would 
seem to be counterintuitive, since it is to be expected that an earthquake should be 
more likely to occur with increasing time in response to an inexorable tectonic loading 
that brings a fault ever closer to its finite threshold of fracture. 

The intuitive interpretation is of course consistent with simple relaxation oscilla- 
tor models of the earthquake process, such as the slip- or time-predictable models. 
But these models should be reconsidered if the stress field is altered on a given fault 
segment due to redistribution derived from earthquakes on nearby fault elements; 
these interactions can cause fluctuations in the stress field, with consequent fluctu- 
ations in the interval times. Knopoff (1996) has proposed that the fluctuations in 
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the interval times between great earthquakes on the San Andreas Fault (Sieh, et al., 
1989) may be associated with stress interactions between the San Andreas Fault and 
other, nearby faults. 

Below, we give a rigorous statistical framework for the derivation of a quantitative 
response to Q. Statistical estimates of recurrence times will be found to be very 
sensitive to assumptions about statistical distributions. Our results confirm, quantify 
and extend the numerical analyses of Davis, et al. (1989) and Ward and Goes (1993) 
and Goes and Ward (1994), by providing an analytic basis for the problem. 



The time to the next earthquake 

Let p(t) be the probability density of the time intervals between earthquakes. If 
the time (now) since the last earthquake is t, what is the probability density function 
P(t') that we must wait an additional time t' until the next earthquake? From 
Bayes' theorem for conditional probabilities, cited in elementary statistics textbooks, 
the probability that an event A, given the knowledge of an event B, is simply the 
quotient of the probability of the event A without constraint and the probability of 
event B: 

P(A\B) = m.. (!) 

Applied to this problem, P(A) = p(t + t') which is the probability that the next 
earthquake will occur at time t' from now, and P(B) = J t °° p(s)ds, which is the 
probability that no earthquake has occurred up to now. Thus 

Pit') - P{t + f) (2) 
F[t) - f t °°p(s)ds' [2) 

which is normalized. 

We calculate the expected time until the next earthquake (t 1 ) as a function of the 
time since the last one. 

(A.) The answer to Q. is given by the sign of ^-r-, if (t) exists. 

From eq. (2), the average expected time to the next earthquake is 

= JTiW™ (3) 

J t °°p(u)du 



By a simple change of variable, 



, = S t °°(u-t)p(u)du ,^ 
!™p{u)du 

We integrate the numerator of (4) by parts and get 

_ IT ds S!° P(u)du 

{} ~ jrp(u)du ■ (5) 

The denominator and numerator of (5) are the familiar first cumulative integral and 
the less familiar second cumulative integral of p(u). For simplicity, we write (5) as 

(0 = -|| (6) 



where f"(t) = pit), i.e. fit) is the second cumulative integral of p(u). Thus 

rlli'\ 

>0 if f(t)f"(t)-[f(t)Y>V . (7a) 



dt 

Equivalent ly, 

^iiH " f /(t)<0 ' (7i) ' 

where we have set f(t) = e~ 9 ^\ The signs are appropriately reversed in the case 
g"(t) > 0. Equation (7b) especially favors an appreciation of the behavior at large 
values of elapsed time t. 

If p(t) is finite at t — 0, we can find yet a third version of (5) which is useful for t 
small. A straightforward expansion for small t shows that 

^ d M =m A-l , (7c) 

where 

A = / ds p(u)du = / sp(s)ds = (t) 
Jo Js Jo 

The second integral on the r.h.s. is obtained by integration by parts; (t) is the average 
(unconditional) time of recurrence between two earthquakes. Let r = ^y, where r 
is the estimate of the waiting time until the next earthquake made immediately after 
the occurrence of the preceding earthquake. We call r the instantaneous estimate of 
t'. Thus, can be rewritten in the simple form 



hm -V-i = l±-i . (7d) 
t^o+ dt r v 1 
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• If the instantaneous estimate r of the waiting time is smaller than the average 
waiting time (t) , the time to the next earthquake increases with increasing time 
since the last one for small t : this reflects the fact that the average waiting time 
(t) is formed by contributions from the distribution over all time, and a value 
of (t) larger than r indicates contributions from the distribution that are larger 
than t at non-zero times; in this case lim t ^ 0+ > 0. 

• If (t) < r, the reverse is true, shorter and shorter time scales are sampled on 
the average as time increases, and the time to the next earthquake decreases 
with increasing time since the last one for small t. 

In particular, if p(0) = 0, then \im t _> 0+ ^§r- = —1, and the time to the next 
earthquake decreases with increasing time the since last one for small t; if however 
p(0) = oo, then lim t _ >0+ ^jr = oo, and the time to the next earthquake increases 
with increasing time since the last one for small t. 

The generalization of (7c) to times other than t = is the criterion (7a), when 
the mean (t) exists. When it does not exist, as for example when the tails of the 
distributions decay slower than t~ 2 , we must compare P(t') as given by (2), with 
p(f). 



Exponential distribution 



In order to develop some intuition, we first consider the exponential distribution, 
which is the familiar case of Poissonian statistics, 

e -t/t 

pit) = -—, t > 0, 
to 

where to is the mean interval time between earthquakes. From eq. (2), 

p-t'/to 

P{t>) = ——. (8) 

Not unexpectedly, the estimate of the time of occurrence of the next earthquake does 
not depend on the elapsed time: the average time from now to the future earthquake 
is to, no matter what the value of t. This case is memoryless; indeed it is the 



6 



only distribution that has no memory. The expected time to the next earthquake is 
(f) = t ; there is no need to invoke the machinery of (7) to derive = 0. The 
Poisson distribution is the unique case g"(t) = 0, which gives the same result. 

The exponential distribution is the fixed point of the transformation p(t) — > P(t'), 
i. e. it is the solution to the functional equation 

P(0=P(0- (9) 

To verify that (8) is the solution to (9), differentiate (2) with respect to t, and 
substitute in (9). We get 

- P (tW) = *<^>. (io) 

Take the Laplace transform of (10) with respect to t', with l/i the transform variable. 
The result (8) follows. Thus the exponential distribution is the fixed point of (2). 

We restate these results: Except for the Poisson distribution, all statistical distri- 
butions must have an average time from now to the future earthquake that depends 
on the time since the last earthquake. If the long-time tail of the function f(t), de- 
fined as the integral of the integral of the distribution p(t), falls off at a rate that 
is faster than exponential, the expected time to the next earthquake is reduced, the 
longer the elapsed time since the last, and vice versa. The Poisson distribution is the 
cross-over between the two states. The exponential case has neither a positive nor a 
negative response to Q., since the time since the last earthquake has no influence on 
the time of the next. 



Other Conditional Distributions 

We calculate the expected time to the next earthquake for several examples of 
statistical distributions p(t) with memory. We illustrate the results in figure 1 by 
displaying the average time to the future earthquake t' plotted against the time since 
the last earthquake t for selected distributions p(t); the positive or negative slopes of 
the curves give the answers to question Q. The details of the calculations are given 
in the Appendix. 

The analytical results are summarized in Table 1. The times in the table are 
scaled by a characteristic time to for a given distribution; the precise definition of 
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t for each distribution is given in the Appendix. In general, to is of the order of 
the mean time between earthquakes, if it is not so exactly. In most cases we can 
give an answer A. that is valid over the entire range of elapsed times since the last 
earthquake. In some cases, we can only give the answer for the limits of short and long 
time since the last earthquake. The table is arranged to favor the limiting responses, 
even though we may have the complete solution. In the case of the exponential 
distribution, the response is neutral as we have already discussed. In the cases of 
the periodic and uniform distributions, the response is only meaningful for times up 
to t ] in the case of the Gaussian distribution, the response is meaningful for long 
elapsed times, as we discuss in the appendix. In the cases of the Weibull distribution 
with m > 1 and the semigaussian distribution, the response can be proved to be 
negative for both short and long elapsed times, and can be inferred to be negative 
for all elapsed times. For the Weibull distribution with m < 1 and the power-law 
distribution, the response is positive for both short and long elapsed times, and can 
be inferred to be positive for all elapsed times. For the lognormal, power-law and 
truncated power-law distributions there is one response for short elapsed times and 
the opposite response for long elapsed times, with the implications of a crossover 
and hence neutral response at an intervening time scaled by to] the lognormal and 
truncated power-law distributions have opposite responses to each other in the long 
and short time regimes. 

Since the truncated and ordinary power-law distributions give opposite results for 
long elapsed times, it follows that the answers to Q. are unstable with respect to the 
presence or absence of a cutoff in the distributions. It is by definition problematical 
that a presumed existence of a cutoff can be identified from a finite set of observations 
of interval times: there is no guarantee that a presumed cutoff will not disappear 
with a future observation of a longer interval between earthquakes. Thus a positive 
response to the question for long elapsed times since the last earthquake is only 
conjectural, i.e. it is only as strong as one knows the distribution to times longer 
than have been observed, which is impossible. Of course, the distribution can always 
be postulated a priori, as in the numerical examples of Davis, et al. (1989) and Ward 
and Goes (1993), but the postulate does not ensure that it represents nature. 
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Estimate of to 



Suppose that we do not know a priori the characteristic time to of time intervals 
between successive earthquakes. We then have to estimate it from a finite suite of 
observations of interearthquake time intervals. Assume that (n — 1) observations 
of time intervals ti, t 2 , ••• , t n -i, are made precisely; we ignore here the additional 
problem of the uncertainties in the time intervals that occur for historical earthquakes; 
this can be treated by standard statistical methods (Sieh, et al., 1989). Suppose that 
the time since the last event is t. Then, in the case of the Poisson distribution 
pit) = e * *° , the standard maximum likelihood method gives the estimate of t as 

En — 1 

the value which maximizes -^e *o i.e. 

1 / \ 

*o = -(*+E*ij- (") 

n \ j=1 / 

Thus, even for the Poissonian case, the use of the information that no event has 
occurred since t, gives an estimate of the average recurrence time to for the next event 
that increases with t ! The Poisson distribution is memoryless only if its parameter 
to is known a priori. 

The calculation (11) can be generalized for the other distributions discussed above. 
Our previous results must thus be reconsidered if the parameters of the distributions 
are themselves not known precisely but are estimated using presently available infor- 
mation. This does not pose any difficulty in principle but must be addressed case 
by case. This simple calculation demonstrates the sensitivity of the "prediction" to 
the assumptions concerning what is really known and what is only inferred from the 
data. 

Summary 

These observations can be summarized as follows: 

• The Poisson or exponential distribution is memoryless and the expected time 
until the next event is independent of previous observations and of the elapsed 
time since the last earthquake. The exponential thus acts as a fixed point in the 
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space of distributions of the transformation (2), and sits at the boundary be- 
tween the positive and negative classes of memory, i.e. at the boundary between 
positive and negative responses to Q. Any statistics of the fluctuations of re- 
currence times that is different from Poissonian entails the explicit assumption 
of a memory. 

• a) Any distribution that falls off at large time intervals at a faster rate than 
an exponential, such as the periodic, quasiperiodic, uniform and semigaussian 
distributions and the Weibull distribution with m > 1, has the property, "the 
longer it has been since the last earthquake, the shorter the expected time until 
the next". The truncated power-law distribution for times close to the cutoff 
time also has this property. 

b) Any distribution that falls off at large time intervals at a slower rate than 
an exponential, such as the Weibull distribution with m < 1, the unbounded 
lognormal and power-law distributions, and the truncated versions of these 
distributions for times remote from the cutoff, have the property, "the longer it 
has been since the last earthquake, the longer the expected time till the next" . 

• a) All distributions that have an instantaneous expectation time interval be- 
tween earthquakes smaller than the average waiting time between earthquakes 
have the property of an increasing time to the next earthquake for an increasing 
time the since the last one, for short times since the last one. This includes the 
cases p(0) = oo. 

b) All distributions that have an instantaneous expectation time interval time 
between earthquakes larger than the average waiting time between earthquakes 
have the property of a decreasing time to the next earthquake for an increasing 
time since the last one, for short times since the last one. This includes the 
cases p(0) = 0. 

• Caution should be exercised in the use of statistics of fluctuations of interval 
times deduced from data sets that describe only the distributions for short time 
intervals between earthquakes. This is because of the strong dependence of our 
result for long times, and in some cases for short times as well, on the properties 
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of the tails of the distributions as well as on the values of the parameters of the 
distributions. 

• The estimate of the time until the next earthquake depends on a precise estimate 
of the tail of p(t) and is unstable with respect to presently available data for 
the recurrence of large earthquakes. Even a finite sampling of the Poisson 
distribution will lead to an estimate of the time to the next earthquake that 
increases with increased time since the last one. 

Thus the positive response of Davis, et al. (1989) to Q., "the longer it has been 
since the last earthquake, the longer the expected time till the next", is shown to 
arise from the use of a distribution that decays slower than an exponential and which 
is unbounded for large time intervals; the result is valid for other distributions as 
well. If a slowly decaying law itself undergoes a transition at even longer intervals to 
a more rapidly decaying law, as in the extreme case of a distribution with a cutoff, 
one can expect that eventually the next earthquake will become more and more 
probable. The response to the question (Q.) is also related, in part, to the finiteness 
of the number of observations of time intervals between earthquakes that give the 
estimate of the distribution p(t); the extrapolation of the estimate of the distribution 
to its asymptote for very large time intervals is exceedingly dangerous, since this 
is procedure is likely to be based on few, if any, observations. The results of this 
exercise suggest that caution be used in the extrapolation of statistics deduced from 
short time-scale data sets to long time-scales. 
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Appendix 



We calculate the expected time to the next earthquake for several examples of 
statistical distributions p(t) with memory. 

A. Periodic distribution 

The simplest of the distributions with memory is the periodic distribution, 

p(t)=8(t -t). 

By inspection, we have 

P(t') = 5(t -t'-t), (t')=t -t. 

Without invoking the generalized machinery, — — 1. In this simplest of cases, 
the expected time of the forthcoming earthquake decreases as the elapsed time since 
the preceding earthquake increases. Extension to quasiperiodic cases can be made. 

B. Uniform distribution 

The uniform distribution is 

p(t) = i < t < t , 

to 

where to is the maximum interval between earthquakes and to/2 is the mean. From 
eq. (2), we get 

P{t') = — !— < t' < t - 1. (Al) 

to — t 

-P(t') is independent of t', i.e. it is itself a uniform distribution, but its value is 
dependent on t. The probability that an earthquake will occur at any time in the 
future up to t increases as the time since the last earthquake increases, and becomes 
infinite as t — > t , which simply expresses the intuitive result that the event will 
occur with certainty before t . It is easy to see that the average time to the future 
earthquake from the present is (t') = |(t — t). The negative value of d{t')/dt gives 
the answer (A.): the expected time to the next earthquake decreases with increasing 
time since the last earthquake. In figure 1, we show the average time to the future 
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earthquake plotted against the time since the last earthquake; both coordinates are 
normalized by the mean time between earthquakes, to/2- The linear relationship 
between (£') and t is strongly curved on the log-log plot. 

In Fig. 2, we display the probability that the next earthquake will occur at 
time t' from now. We show the unconditional probability, i.e., the probability as 
though we knew the distribution of intervals p(t) but did not know the time of the 
last earthquake. We also show the (conditional) probability of an earthquake in the 
future knowing that the last earthquake took place at time 1.33to in the past. In the 
latter case, no earthquake can occur after 0.67t from now; according to (Al), the 
probability of occurrence of the future earthquake is higher by a factor of 3 than the 
unconditional probability, and is independent of the time of the future earthquake. 

C. (Semi)-Gaussian distribution 

The Gaussian distribution is 



where t is both the mean and the most probable time interval of earthquake recur- 
rence; the standard deviation is a. The Gaussian distribution has a finite probability 
that the next earthquake will occur before the preceding one. The drawback is minor 
if a <C t , a condition that describes a nearly periodic distribution; we have consid- 
ered the periodic case above. To restrict the problem to cases of positive t, we could 
truncate the distribution at t — 0, however this leads to messy mathematics; for large 
times, the drawback is minor. 

For the simpler problem t ^> t , we use the approximation 



which is the semi-gaussian distribution, i.e. it is a Gaussian centered at t — 0, the 
mean time interval of earthquake recurrence is to an d the most probable time for 
recurrence is of course zero for this distribution. Formally, eq. (2) yields 
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(t-tpf 
e 2 CT 2 



pit) 




t > 




(A2) 
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where erfc(x) is the usual complementary error function. 

If the elapsed time since the last earthquake is very large, t ^> to, we can use 

the first term of the asymptotic expansion of the second cumulative integral which is 

-A 

fit) ~ ^i? 1 - Then g(t) ~ t 2 + O(logt), g"{t) > for large t, and from (7b) 



t 2 

d(t') 
dt 



< 



For short times, we use (7c) with p(0) = ^ and A = t and get 

- 1 = -0.363. 



TVto 

d(t'} ^ 2 



dt 7T 

Thus in both the short and long time limits, the expected time until the next earth- 
quake decreases as the time since the last one increases. 

D. Lognormal distribution 

The lognormal distribution is 

m = = 'e-^i^. (.43) 

which is similar in shape to the Rayleigh distribution (see below) near t — but has 
a much more slowly decaying tail for large times. In this case, t is the median time; 
the mean time is (t) = t e~; the most probable time is to e ~ CT • From (A3), we can 
write 

^W = (log(f) + ^ 2 ) 2 + 0(loglog(f)), 

whence 

/W-|(l-^ 2 -log^). 

For t > t , g" < and hence > 0. For t < t , g" > and < 0; 
alternatively, we note that p(0) = and hence from (7c), ^r- < 0, which is the same 
result. Thus the lognormal distribution has a crossover in response to Q. 

We express these results graphically. For the case a = t we display P(t') for times 
t = 2to and t = 5t (Fig. 3) together with the unconditional log-normal distribution 
p(t). From figure 3a, we see that P(t') is significantly smaller than p(t') for times 
comparable to the elapsed time t but P(t') is, as expected, larger than p(t') at large 
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times (see extension of Fig. 3a to long times in Fig. 3b). This is a small effect 
for t = 2t but is much stronger for t = 5t and all the more so if t increases even 
more. Thus, numerically as well as analytically, for early elapsed times t that are 
comparable to the peak of the distribution, the longer the elapsed time since the last 
event, the shorter the time until the next event; but for large elapsed times since 
the last earthquake, the longer the time since the last event, the longer the time 
until the next one. In the lognormal case, it is correct to state that "the longer it 
has been since the last earthquake, the longer the expected time until the next" but 
only for elapsed times greater than times of the order of the characteristic time. The 
lognormal distribution is an example of a case that has one answer to Q. for short 
times since the last earthquake and the opposite answer for long times. Of course, 
the crossover takes place at elapsed times that are of the order of the characteristic 
time to- Note that the probabilities for long times in the future are, as might be 
expected, very small. The log-normal case is an example of a distribution with a tail 
that decays slower than an exponential. 

To illustrate more concretely the properties of a system that has a crossover 
response to Q., we concoct the distribution 



which probably has no redeeming virtue in nature, but has the property that it has 
easily calculable integrals. It is evident that this distribution has both a long-time 
tail that decays slower than exponential, and the property p(0) = 0, as in the case 
of lognormal distribution. Thus we are guaranteed that there is a crossover in the 
response to Q. More precisely, the criterion function (7a) is f(t)f"(t) — [f'(t)] 2 = 
j(x 3 + Ax 2 + Ax — A)e~ 2x where x = y^. The criterion has the crossover in sign 
between x small and large at ^ = 0.3532. The lognormal distribution has similar 
properties. 
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E. Weibull distribution 

The Weibull distribution is 

p(t) = mto m t m - x e~^ ] , < t < oo, m > 

having a most probable value (m — l)™t , a mean a(m)t where a(m) = J °° e~ tm dt, 
and median (log2) 1 / m to- Values of m smaller than 1 correspond to p(t) decaying 
slower than an exponential for large t and give the so-called stretched exponential 
distributions, while values of m larger than 1 lead to a decay that is faster than 
exponential. For large m, pit) approaches a delta function centered on t , i.e. to the 
periodic distribution we have considered above. Ward and Goes (1993) and Goes and 
Ward (1994) have studied the degree of earthquake clustering as a function of m; in 
their notation, v — 1/m. Equation (2) yields 



f f t) „ e ° = e -(.i- r-^-i)io 9 i- 



P{t') = mt m (t + tT^e ^ . (AA) 

The first term of the asymptotic series for /(£), which is the second cumulative integral 
of p(t), is 

-(— y 

(±)3(m-l) 

Evidently g(t) ~ t m + 0{logt) and hence g"(t) ~ m{m - l)t {m ' 2 \ Thus if < m < 1, 
then g"(t) < 0, and ^ > for large t, while if m > 1, then g"(t) > 0, and ^ < 
for large t. For small t, p(0) = 0,m > 1 and it follows from (7c) that ^ < 0. For 
cases m < 1, p(0) = oo, and from (7c), > 0. Thus there is a reversal in response 
between the cases m < 1 and m > 1, in agreement with the result of Ward and Goes. 

In figure 4a, we exhibit the interesting subcase of the Weibull distribution with 
m = 2 which is appropriate for rectified Gaussian noise and is known as the Rayleigh 
distribution. The Rayleigh distribution has a tail with similar properties to that 
of the Gaussian, and decays faster than an exponential. We show Pit 1 ) for times 
t — t and t = 2t . It is clear that P(t') has a progressively shrinking width to the 
origin as t increases, i.e. the expected time decreases as the waiting time t increases. 
The answer A. is negative for both short times and long times, the latter property 
evidently connected with the rapid fall-off in p(t) for large t. 
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The opposite situation is found in the case m < 1; in figure 4b we plot the case 
m = 1/2 and show P(t') for times t = to, t = 2t and t = 10t , as well as p(t). It 
is clear that P(t') lies well above p(t') at long times t' > to, and all the more so as 
t increases. Thus the longer we wait, the longer the time to the next event, in this 
case. 

F. Power-law distribution 

The unconditional power-law distribution is 

p(t) = 0, < t < t , (-45) 

aft \ 

p{t) = I — ) t <t< OO, 0</!<OO 

to \to/ 

The characteristic time scale to is proportional to the mean -^t for p > 1 and the 
median t 2^; for p < 1, the mean is infinite. This is an example of a distribution 
with a waiting time and has been used in the case p = 1/2 in short-term earthquake 
prediction calculations by Kagan and Knopoff (1981, 1987). 

For p > 1, we evaluate for this case by applying the criterion function (7a), 
which gives (jff 1 * Thus ^ > for all p > 1 and all t > t . 

For p < 1, (t) diverges and the criterion (7a) cannot be used. In this case, we 
have to examine the conditional distribution (2) directly, and compare it with the 
unconditional distribution. This general method is also applicable to the case p > 1. 
Substitution of (A5) in eq. (2) yields 

Pit') = + iA6) 

t t 

Formula (A6) is almost the same as the unconditional distribution, except for the 
additional 1 in the parentheses; the two expressions are identical in the limit t' ^> t 
except for obvious scaling factors. Thus in this limit, the distribution P(t') is also 
a power law with a characteristic scale given by the waiting time t, instead of t for 
the unconditioned pit'). Thus the longer it has been since the last earthquake, the 
longer the expected time until the next, for all cases p > 1. 

Figure 5a shows Pit') for t = 10t and t = 100t together with pit) for an 
exponent p = 3 (for this choice, p(t) possesses a finite mean and variance). We 
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observe the asymptotic power-law behavior of Pit') at times t' > t with amplitude 
significantly larger than p(t'), showing the enhanced probability for large conditional 
waiting times. Figure 5b shows P(t') for t = 10t and t = 100t together with pit) 
for the threshold case of exponent /i = 1; in this case, the mean and variance are not 
defined. This is an illustration of a power law with a very weighty tail. The behavior 
of Pit') is qualitatively similar to the previous case. 

G. Truncated power-law distribution 

Except for the uniform and periodic distributions, we have considered thus far 
only distributions of fluctuations in interval times that extend to infinity. In these 
cases of distributions with long-time tails, there is a finite but small probability that 
a second earthquake will occur after a very long time interval after the first. If the 
distributions describe the seismicity of a region, rather than that of an individual 
fault, the very long time intervals imply very large accumulations of deformational 
energy and hence very large fracture sizes. To avoid the problems of earthquake sizes 
greater than the size of a given region, we consider a cutoff in the distributions pit) 
(Knopoff, 1996). To demonstrate the influence of a cutoff, we restrict the previous 
case to 

P(t) = ' to <t< t max 

= < t < t 

which is normalized. Substitution in eq. (2) yields 

m = JtTW^ t-'-iZL ' {A7) 

For t <C t max , the second factor of (A7) is t^, which is the same as letting t max — > oo. 
Thus we recover the previous case of the simple power law without truncation. 

The interesting regime is found when t is not very small compared to t max . Con- 
sider the case t — > t max . From eq. (A7) we get 

Pit') * , (A8) 

tmax £ 

which is independent of t' and becomes very large as t — > t max . This case is identical 
to that of the uniform distribution (Al) above. Thus it is not unexpected that the 
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longer we wait, the shorter will be the expected time until the next event. Without 
truncation, the result is reversed. There is a crossover between the truncated and 
untruncated cases as illustrated in figure 6. In the figure, we take \x = 3 as in figure 
5a, t max = 100£ and show P(t') for t = 10t , 90t and 98t - For t = 10t , P(t') is 
found to be much larger than p(t') in the tail, as in the previous untruncated case. 
For t = 90*o, P(t') is defined only for < t' < I0t . In agreement with (A7), we see 
that P(f) becomes almost constant and close to For t = 98t , P(t') is defined 
only for < t' < 2t and is close to This illustrates the crossover from a longer 
expected time when t is small, to a shorter expected time as t approaches t max - 

Since the truncated power-law distribution is very close to a uniform distribution 
for times near t max , we expect that the truncated lognormal distribution and trun- 
cated Weibull distributions with /i < 1 will also have a shorter time until the next 
event, the longer we have been waiting for an earthquake to happen; thus we expect 
a crossover in the response to Q. between these truncated and untruncated cases as 
well. 
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TABLE 1 : RESPONSE TO Q. 



Distribution 


Short times 


Long Times 


Exponential 








Periodic 


- 




Uniform 






Gaussian 






Semigaussian 






Lognormal 




+ 






+ 


Weibull (m > 1) 






Rayleigh (m=2) 






Weibull (m < 1) 


+ 


+ 


Power Law 




+ 


Truncated Power Law 


+ 
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Figure Captions 



Figure 1. Expected time of the next earthquake as a function of the elapsed time 
since the last one. The response to Q. is given by the position of the curve with 
respect to 1. 

Figure 2. Uniform distribution: p(t') and P(t') for t = (4/3)t - 

Figure 3. Lognormal distribution with a = t$: p(t') and P{t') for t = 2t and 
t = 5t . a) short times; b) long times. 

Figure 4. Weibull distribution. 

a) m = 2 (Rayleigh distribution), corresponding to a tail decaying faster than an 
exponential. P[t') is shown for t — to and t = 2t together with pit'). 

b) m — 1/2 corresponding to a tail decaying slower than an exponential. P(f) is 
shown for t = t , t = 2t and t = 10t , together with p(t'). 

Figure 5. Power- law distribution. 

a) /i = 3; p(t) possesses a finite mean and variance. P(f) is shown for t = 10to 
and t = 100t together with p(t'). 

b) /i — 1; the mean and variance are not defined. P(t') is shown for t = 10t and 
t = 100t together with p(t'). 

Figure 6. Truncated power-law distribution with /i — 3 and t max = 100t - Pit') 
is shown for t = 10t , 90t and 98t - 
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